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We report on the measurement of the Cabibbo-Kobayashi-Maskawa angle 7 through a Dalitz 
analysis of neutral D decays to K s tv~it + in the processes 

D* — » _D7r°,Z)7, with the BABAR detector at the SLAC PEP-II e + e~ asymmetric-energy collider. 
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I. INTRODUCTION AND OVERVIEW 

The angle 7 of the unitarity triangle is the phase of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix [l[ de- 
fined as 7 = arg [— V ud V* b /V cd V* b ], which corresponds to 



Kte- 



the phase of the element V* b , i.e. V ub 
the Wolfenstein parameterization The precise mea- 
surement of the angle 7 is a crucial goal of the physics 
program at the B-factories, however, it is also one of the 
most difficult to achieve. 

Among all methods proposed to extract 7, only those 
using B T — > D°K T decays are theoretically clean be- 
cause the main contributions to the amplitudes come 
from tree-level diagrams (the symbol D° indicates ei- 
ther a D° or a D° meson) . The interference between the 
color allowed B~_— > D°K~ (b — > us) and the color sup- 
pressed B~~ -^_D°K~ (b — » ucs) transitions |3j, when 
the D° and D° are reconstructed in a common final 
state d, [H, H, 0], introduces a relative phase 7 in the 
decay amplitude. The sensitivity to 7 depends on the 
magnitude of the ratio of the b — > ucs amplitude with 
respect to the 6 — > ys one, re, which plays a key role on 
the ability to measure 7 at the B-factories. Theoretical 
expectations, consistent with current experimental lim- 
its, give r B »| V^^/V^VS I c F ~ 0.1, where c F ~ 0.2 
is the color suppression factor. 

When the D° is reconstructed in a 3-body final state 
like K®tt~tt + , the interference between doubly-Cabibbo 
suppressed, Cabibbo allowed and CP-eigenstate ampli- 
tudesprovides strong phases to ensure the sensitivity to 
7 0, 0] • The angle 7 can then be extracted through an 
analysis of the distribution of the events in the D° Dalitz 
plane. 

Assuming negligible effects from D° — D° mixing [8| 
and CP asymmetries [9] in D decays, the B T — » 
D(*)°K^, with D*° -> D°tt q ,D°j, D° -► K°it-it + decay 



chain rate r^\m?_ 



i 2 ,) can be written as 



T%\m 2 _,ml) « I^tI 2 + 4* )2 |^ ± | 2 + 

2e {4*' Re[A DT A* D± ] + lm[A DT A D± }} , (1) 
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where m 2 _ and m 2 ^ are the squared invariant masses 
of the K ( *tt~ and K®tt + combinations, respectively, 
_4d=f = ^n(m^,m|), with Ad- (Ad+) the ampli- 
tude of the D' 







K° s 7r-n+ (D° -> K^n+n~) de- 
cay. We introduce the CP ( cartesian ) parameters (Tpj 



„(*) _ J,*) 



cos(5 B 



(*) 



=F 7) and y 



(*) 

verifying 



(*) 



(*) _ 

Here. 



„(*) . 



b sin((5 B 



(*) 



„(*) 



T 7), 
is the magni- 
tude of the ratio of the amplitudes A(B~ — > D(*^°K~) 
and A(-B~ — > Z>W°if _ ) and 5^ is their relative strong 
phase. The factor e in Eq. (JTJ) takes the value —1 for the 
decay B^ -> D*°[D°i\K^ and +1 for all the rest. This 
relative sign arises due to parity and angular momentum 
conservation in the Z)M° decay, and the different CP 
content of D°7 with respect to D°tt° [111 ]. 

Equation |TJ also applies to B T — ► D°K* T decays, 
with the replacements r B - 

nr s cos(£ s =F 7)7 and -> y ST 
ing + j/g-p = n 2 r 2 . Here, the parameter k accounts 
for interference between resonant and non-resonant K* 
decays, as a consequence of the natural width of the K*, 
with < k < 1 [12]. This general parameterization also 
accounts for variations of r s and S s within the K* mass 
window, and for efficiency variations as a function of the 
kinematics of the B decay. 



8g ► S Sl X~ > Xgzf: — 

nr s sin(<5 s Tl), verify- 



II. DATA SAMPLE AND EVENT SELECTION 

The analysis for B~ -► D^°K~ {B~ -> D°K*~) de- 
cays [|| is based on a sample of approximately 347 (227) 
million BB pairs collected by the BABAR detector [13| at 
the SLAC PEP-II e + e~ asymmetric-energy storage ring. 
For each signal B decay channel we also reconstruct its 
own control sample, B~ — > D^°tt^ (B~ — > D°a^). 

The reconstruction and selection criteria are described 
in detail elsewhere [13, E3, [3- B meson candidates 
are characterized by using the energy difference AE, the 
beam-energy substituted mass jties, and a Fisher dis- 
criminant T to separate e + e~ — > qq, q = u, d, s, c (contin- 
uum) and BB events 0. If both S~ -> l)* !!? ^ ]^- 
and B~ — * _D*°[£)°7]iir _ candidates are selected in the 
same event, only the B~ -> L>* [I? 7r ]A^ is kept. The 
cross-feed among the different samples is negligible ex- 
cept for B~ — > D*°[D°7]i ; r _ , where the background from 
B- -» £)* [I) 7r ]A:- is about 5% of the signal yield. 
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This contamination has a negligible effect on the CP pa- 
rameters. 

The reconstruction efficiencies are 15%, 7%, 9%, and 
11%, for the B~ -> D°K~, B~ -> ^[DV^-, 
S~ -> £>*°[D°7]ir-, and B~ -> D°K*~ decay modes, 
respectively. Figure [T] shows the mES distributions after 
all selection criteria, for \AE\ < 30(25) MeV, for B~ -> 
£)(*)" R-(B- -> D°K* ). The largest background con- 
tribution comes from continuum events or BB decays 
where a fake or true D° is combined with a random 
track. Another source of background for B~ -> D^°K~ 
is given by — > D (*)°7r~ decays where the prompt pion 
is misidentified as kaon. These decays are separated from 
the signal using their different AE distribution. 




m ES (GeV/c 2 ) 



FIG. 1: Distributions of m E s for (a) B~ -> D°K~, (b) 
B~ -> 5*°[5 7r°]if _ , (c) S- -> D* ^ ^", and (c) 
_B~ — > D°K*~ . The curves superimposed represent the over- 
all fit projections (solid black lines), the continuum contribu- 
tion (dotted red lines), and the sum of all background com- 
ponents (dashed blue lines). 



III. THE D° -> A'Svr 7T+ DECAY MODEL 



The D° 



7r + decay amplitude _4o(m? 



-) 



is determined from an unbinncd maximum-likelihood 
fit to the Dalitz plot distribution of a high-purity 
(97.7%) tagged D° sample from 390328 D*+ -> D°tt+ 
decays reconstructed in 270 fb _1 of data, shown in 
Fig. [21 Our phenomenological reference model to de- 



scribe Ad(w: 



uses a sum of two-body amplitudes 



(subscript r) and a non-resonant (subscript NR) contri- 
bution, 

A D (mL,m%) = ^2a r e i<t,r A(m 2 _,ml) + a NR e^ NR , (2) 



where the parameters a r (onr) and <p r (</>nr) are the 
magnitude and phase of the amplitude for component r 
(NR). The function A r = F r x T r x W r is the Lorentz- 
invariant expression that describes the dynamic proper- 
ties of the D° meson decaying into K®tt~tt + through an 
intermediate resonance r, as a function of position in the 
Dalitz plane. Here, F r is the Blatt-Weisskopf centrifugal 
barrier factor for the resonance decay vertex (l6j with ra- 
dius R — 1.5 GeV -1 (0.3 fm), T r is the resonance prop- 
agator, and W r describes the angular distribution in the 
decay. For T r we use a relativistic Breit-Wigner (BW) pa- 
rameterization, except for r — p(770) and p(1450) where 
we use the functional form suggested in Ref. [17] . The 
angular dependence W r is described with the helicity for- 
malism as shown in [18I| 1 . Mass and width values are 
taken from [ljj], with the exception of ^Cq(1430) + taken 
from [2(|. The model consists of a total of 13 resonances 
leading to 16 two-body decay amplitudes and phases (see 
Table[T]), and accounts for efficiency variations across the 
Dalitz plane and the small background contribution. All 
the resonances considered in this model are well estab- 
lished except for the two scalar 7T7t resonances, a and 
cr', whose masses and widths are obtained from our sam- 
ple [2l| . Their addition to the model is motivated by an 
improvement in the description of the data. 

The possible absence of the a and a' resonances is con- 
sidered in the evaluation of the systematic errors through 
the use of a K-matrix formalism to parameterize the 
7T7T S-wave states. The K-matrix method provides a di- 
rect way of imposing the unitarity constraint of the scat- 
tering matrix that is not guaranteed in the case of the 
BW model and is suited to the study of broad and over- 
lapping resonances in multi-channel decays, avoiding the 
need to introduce the two a scalars, 



A D (m 2 



-)=F 1 (s) 



"' ' 

r^nir S=0 



l ^A r {r 



-), (3) 



where Fi(s) — J2j [I ~ *^( s )p( s )]i/ Pj( s ) is the contri- 
bution of tttt S-wave states. Here, s = m 2 _ , , I is the 
identity matrix, K is the matrix describing the S-wave 
scattering process, p is the phase-space matrix, and P is 
the initial production vector [22| . The index j represents 
the j th channel (1 = tttt, 2 = KK, 3 = multi-meson (23|, 
4 = rj-q, 5 = T]7]'). The K-matrix parameters are obtained 
from a global fit to the available tttt scattering data below 
1900 MeV/c 2 [13], while the initial production vector is 
obtained from our fit to the tagged D — » K®tt~tt + data. 



IV. CP FIT RESULTS AND SYSTEMATIC 
UNCERTAINTIES 

Once the decay amplitude ^(m?.,™^) is known it 
can be fed into Eq. |TJ). The extraction of the CP- 



1 The label A and B should be swapped in Eq. (6) of [IS 
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TABLE I: Complex amplitudes a r e l ^ T and fit fractions of the 
different components (Kstt~ , Kstt + , and ir + ir~ resonances) 
obtained from the fit of the D° — > Ksn^TY + Dalitz distribu- 
tion from D* + — > D°ir + events. Errors are statistical only. 
The fit fraction is defined for the resonance terms as the in- 
tegral of a%\A r {m 2 _, m+)| 2 over the Dalitz plane divided by 
the integral of \Ar>(rn^., m+)| 2 . The sum of fit fractions is 
119.5%. A value different from 100% is a consequence of the 
interference among the amplitudes. 



Component 


Re{a r e i ' t ' r } 


Im{a r e i4,r } 


Fraction (%) 


K* (892)~ 


-1.223 ±0.011 


1.346 ±0.010 


58.1 


tf *(1430)- 


-1.698 ±0.022 


-0.576 ± 0.024 


6.7 


(1430)" 


-0.834 ±0.021 


0.931 ±0.022 


6.3 


iT(1410)- 


-0.25 ±0.04 


-0.11 ±0.03 


0.1 


^*(1680)" 


-1.285 ±0.014 


0.205 ±0.013 


0.6 


7y*(892)+ 


0.100 ±0.004 


-0.127 ±0.003 


0.5 


K *(1430)+ 


-0.027 ±0.016 


-0.076 ±0.017 


0.0 


^ 2 *(1430)+ 


0.019 ±0.017 


.177 ±0.018 


0.1 


p(770) 


1 





21.6 


^(782) 


-0.0219 ±0.0010 


0.0394 ± 0.0007 


0.7 


M1270) 


-0.699 ±0.018 


0.387 ±0.018 


2.1 


p(1450) 


0.25 ± 0.04 


0.04 ±0.06 


0.1 


Non-resonant 


-0.99 ±0.19 


3.82 ±0.13 


8.5 


/o(980) 


0.447 ± 0.006 


0.257 ± 0.008 


6.4 


/o(1370) 


0.95 ±0.11 


-1.619 ±0.011 


2.0 


a 


1.28 ±0.02 


0.273 ± 0.024 


7.6 


a 


0.290 ±0.010 


-0.066 ±0.010 


0.9 



violating parameters x^* and (CP fit) is then per- 
formed through a simultaneous maximum likelihood fit 
to the r_(rn?_,m\) and 1+ m\) Dalitz plot distri- 

butions for B~ -> £>W°K" _ and B+ -> D^ Q K+ decays, 
respectively. A similar fit is performed for B T — > D°K* T 
decays. Different background components are consid- 
ered: continuum, B~ -> D^°tt- for B~ -> D^°K~, 
and BB. The likelihood function uses the Dalitz plot dis- 
tribution (after correction for efficiency variations), TOeSi 
AE, and J-, with shapes determined directly from the 
signal and control samples, from both signal and side- 
band regions. Only the shapes for BB background events 
are determined from Monte Carlo simulation. Events 
falling into the continuum and BB background compo- 
nents are themselves divided into events with a real or 
a fake (combinatorial) D°. We finally account for the 
correlation between the flavor of true D° mesons and the 
charge of combinational charged kaon. 

We find 398 ± 23, 97 ± 13, 93 ± 12, and 42 ± 8 sig- 
nal events, for B~ -> D°R-, B~ -> D*°[D°tt q }K- , 
B~ -> D*°[D°j}K~ , and B~ -> D°K*~, respectively, in 
agreement with expectations based on measured branch- 
ing fractions and efficiencies estimated from Monte Carlo 
simulation. The results for the CP-violating parameters 

x^\ t/£* , x s= p, and y s zp, are summarized in Table ITT1 The 
only non-zero statistical correlations involving the CP 
parameters are for the pairs (x_,y_), (x + ,y + ), (x*_,y*_), 
(a;*, 2/+), (x s _,y s _), (x s+ ,y s+ ), which amount to -1%, 




m? (GeV 2 /c 4 ) m|„. (GeV 2 /c 4 ) 

FIG. 2: (a) The 5° -> K° s -r--k + Dalitz distribution from 
D*~ — » D°ty~ events, and projections on (b) m\ — m^ < o 7T +, 

(c) m 2 . = m\ ^, and (d) m^+„._. D° — > K° s -k + -k~ from 
D *+ _^ £)%+ events are also included. The curves are the 
reference model fit projections. 



1%, -17%, -14%, -10%, and 2%, respectively. Figured 
shows the one- and two-standard deviation confidence- 
level contours (including statistical and systematic un- 
certainties) in the (x, y) plane for all the reconstructed 
modes, and separately for B~ and B + decays. The sep- 
aration of the B~ and B + contours is an indication of 
direct CP violation. 

The largest single contribution to the systematic un- 
certainties on the CP parameters comes from the choice 
of the Dalitz model used to describe the D° — > K®ir~ir + 
decay amplitude. We use a set of alternative models 
where some resonances are removed or the parameteri- 
zation of the different amplitudes are changed. For the 
7T7T S-wave we use the K-matrix approach described in 
Sec. IIII1 while for the P-wave we change the mass and 
width describing the p(770) within their quoted uncer- 
tainty fl9j |. The uncertainty on the description of the 
Kn S-wave is estimated by floating in our flavor tagged 
D° sample the mass and width of the BW describing 
the if* (1430), and using an additional parameterization 
taken from Ref. [25j with parameters extracted from our 
fit. Since the Kit P-wave is dominated by the if* (892) 
in both Cabibbo allowed and doubly Cabibbo suppressed 
amplitude, the mass and the width of this resonance, 
taken from [19j | in the reference model, are changed to 
the values obtained from our fit to the tagged D° sample. 
The resulting values are consistent with what is found 
in B — > J/^Ktt decays selected in BABAR data. For 
the 7T7T and Kit D-waves, described by the /2(1270) and 
if* (1430) resonances, respectively, we use as alternative 
the formalism derived from Zemach tensors [26| . The dif- 
ference is very small for P- waves but is larger for D-waves. 



4 



TABLE II: CP- violating parameters , y-C , x s ^, and j/ s =p, as obtained from the CP fit. The first error is statistical, the 
second is experimental systematic uncertainty and the third is the systematic uncertainty associated with the Dalitz model. 

CP parameter B~ -» D a K~ B~ -> b*°K~ B~ -> D° K*~ 

x-/xt/x s - 0.041 ± 0.059 ± 0.018 ± O.Olf -0.106 ± 0.091 ± 0.020 ± 0.009 -0.20 ± 0.20 ± 0.11 ± 0.03 

yJ/y*J/y s - 0.056 ± 0.071 ± 0.007 ± 0.023 -0.019 ± 0.096 ± 0.022 ± 0.016 0.26 ± 0.30 ± 0.16 ± 0.03 

x+/x* + /x s+ -0.072 ± 0.056 ± 0.014 ± 0.029 0.084 ± 0.088 ± 0.015 ± 0.018 -0.07 ± 0.23 ± 0.13 ± 0.03 

y+/y* + /y a+ -0.033 ± 0.066 ± 0.007 ± 0.018 0.096 ± 0.111 ± 0.032 ± 0.017 -0.01 ± 0.32 ± 0.18 ± 0.05 
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FIG. 3: Contours at 39.3% (dark) and 86.5% (light) confi- 
dence level (corresponding to two-dimensional one- and two- 
standard deviation regions), including statistical and system- 
atic uncertainties, for (a) (x^,y^), (b) (x^,y^), and (c) 
[x a ^,y B -^.) parameters, for B~ (thick and solid lines) and B + 
(thin and dotted lines) decays. 



Other alternative models are built by removing the Blatt- 
Weisskopf penetration factors (T^ |. removing resonances 
with small fit fractions --KJ(1430), AT* (1680), AT* (1410) 
and p(1450)-, and replacing the running width of the 
BW by a fixed value. As total systematic uncertainty 
associated with the Dalitz model, given in Table [Til wc 
consider the sum square of contributions from each alter- 
native model, where each contribution is evaluated from 
the difference between the CP fit parameters using the 
alternative and the reference models. The dominant con- 
tributions to the overall Dalitz model uncertainty arise 
from the 7T7t and Kit S-waves, and the fixed BW width. 

Experimental systematic uncertainties arise from sev- 
eral sources and can be found in Table lnTl All of them are 
small compared with the statistical precision, and their 
sum is similar to the Dalitz model uncertainty. Other 
possible sources of experimental systematic uncertainty 



are found to be negligible. 



V. INTERPRETATION AND CONCLUSIONS 

A frequentist (Neyman) procedure [l(| [H| has been 
adopted to interpret the measurement of the CP param- 
eters (x^C\y^P) reported in Table UTI in terms of con- 
fidence regions on (7, t\b, 8b, r B , 8 b ). No attempt has 
been done yet to include in the combination the CP pa- 
rameters for B T — > DK* T decays, (x s zf, y s ^)- Figured] 
shows the two-dimensional projections onto the (re, 7) 
and (7-3,7) planes of the one- and two-standard devia- 
tion regions, including statistical and systematic uncer- 
tainties. The figure reveals the two- fold ambiguity of this 



(7,4*i) 



(7 + 180°, 4*1 + 180°), as expected 



method, 

from Eq. |T]), From the one-dimensional projections we 
obtain for the weak phase 7 = (92±41±11±12)°, and for 
the strong phase differences S B = (118 ± 63 ± 19 ± 36)° 
and 8* B = (-62 ± 59 ± 18 ± 10)°. No constraints on 
the phases are achieved at two-standard deviation level 
and beyond. Similarly, for the magnitude of the ratio of 
decay amplitudes r B and r* B we obtain the one- (two-) 
standard deviation constraints r B < 0.140 (r B < 0.195) 
and 0.017 < r* B < 0.203 (r* B < 0.279). No constraint on 
7 is obtained from B T — ► DK* T decays alone, for which 
nr s < 50(0.75) at one- (two-) standard deviation level. 
All these results are obtained considering the statistical 
correlations discussed in Sec. IIV1 while the experimen- 
tal and Dalitz model systematic uncertainties are taken 
uncorrelated. We have verified that accounting for exper- 
imental systematic correlations within a given B decay 
channel, (x T ,y T ), (x^,y^), or (x ST , y ST ), or assuming 
the experimental and Dalitz model systematic uncertain- 
ties between (x Tl y T ) and {x^,y^) fully correlated, has 
a negligible effect on the results. 

In conclusion, BABAR has reached a good precision in 
the measurement of the CP parameters (a;V ,y±. ) but 
the improvement of the statistical and systematic un- 
certainties on 7 also depends on the value of the r B 
parameter (the former scales as l/r B ^). Our last ex- 

pcrimcntal results for r B are somewhat smaller than in 
our previously published analysis |l0j. Therefore, an im- 



proved precision in the determination of 



„(*) 



is funda- 



mental to better constraint 7. In this respect, the anal- 
ysis of more data being recorded by the detector, the 
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TABLE III: Summary of the main contributions to the experimental systematic error on the CP parameters. 



Source 


X- 


V- 


x+ 


y+ 


x*_ 


V- 


x*+ 


y+ 


X s - 


Va- 


x s + 




m,ES, AE, T shapes 


0.002 


0.004 


0.003 


0.004 


0.011 


0.012 


0.008 


0.008 


0.08 


0.12 


0.10 


0.12 


Real D fractions 


0.002 


0.000 


0.000 


0.000 


0.002 


0.003 


0.002 


0.016 


0.03 


0.03 


0.03 


0.04 


Charge- _D° flavor correlation 


0.008 


0.002 


0.002 


0.002 


0.005 


0.005 


0.001 


0.022 


0.03 


0.04 


0.03 


0.05 


Efficiency in the Dalitz plot 


0.014 


0.000 


0.013 


0.001 


0.001 


0.002 


0.000 


0.001 


0.06 


0.04 


0.07 


0.09 


Background Dalitz shape 


0.006 


0.003 


0.001 


0.004 


0.012 


0.015 


0.009 


0.009 


0.04 


0.09 


0.04 


0.09 


Dalitz amplitudes and phases 


0.004 


0.004 


0.004 


0.004 


0.008 


0.008 


0.008 


0.008 


0.01 


0.01 


0.01 


0.01 


B~ -> D*°K~ cross-feed 


0.000 


0.000 


0.000 


0.000 


0.004 


0.001 


0.004 


0.004 










CP violation in D-k and BB bkg 


0.000 


0.000 


0.000 


0.000 


0.005 


0.002 


0.002 


0.005 


0.00 


0.00 


0.00 


0.00 


Total experimental 


0.018 


0.007 


0.014 


0.007 


0.020 


0.022 


0.015 


0.032 


0.11 


0.16 


0.13 


0.18 




FIG. 4: Projections onto the (a) (rs, 7) and (b) (r* B , 7) planes 
of the 3.7% (dark) and 45.1% (light) five-dimensional confi- 
dence level regions, corresponding to one- and two-standard 
deviation intervals, respectively, including statistical and sys- 
tematic uncertainties. 



addition of different B (e.g. B T -> DK*T) an d D (e.g. 
D° —> 7r 7r~7r + , KgK ~K + ) decay channels, and the com- 
bination with other methods 0, H will be helpful. As- 
suming tb = 0.1 it will be possible to measure 7 with 
~ 10° error with a 1 ab _1 data sample, which is within 
the reach of the BABAR experiment. 
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